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THE 7-RAY EMISSION REGION IN THE FRJI RADIO GALAXY 3C 111 
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ABSTRACT 

The Broad Line Radio Galaxy 3C 111, characterized by a Fanaroff- Riley II (FRII) radio morphology, 
is one of the sources of the Misaligned Active Galactic Nuclei sample, consisting of Radio Galaxies 
and Steep Spectrum Radio Quasars, recently detected by the -Fermi-Large Area Telescope. 

Our analysis of the 24-month 7-ray light curve shows that 3C 111 was only occasionally detected 
at high energies. It was bright at the end of 2008 and faint, below the Fermi-Large Area Telescope 
sensitivity threshold, for the rest of the time. A multifrequency campaign of 3C 111, ongoing in the 
same period, revealed an increase of the mm, optical and X-ray fluxes in 2008 September-November, 
interpreted by Chatterjee et al. (2011) as due to the passage of a superluminal knot through the jet 
core. The temporal coincidence of the mm-optical-X-ray outburst with the GeV activity suggests a 
co-spatiality of the events, allowing, for the first time, the localization of the 7-ray dissipative zone 
in a FRII jet. We argue that the GeV photons of 3C 111 are produced in a compact region confined 
within 0.1 pc and at a distance of about 0.3 pc from the black hole. 

Subject headings: galaxies: active galaxies: jets gamma-rays: galaxies: individual (3C111) 



1. INTRODUCTION 

The Large Ar ea Telescope (LAT) on board the Fermi 
7-ray satellite (|Atwood et al. 20 09) has recently de- 
tected a handful of Misaligned Active Galactic Nuclei 
(MAGN), i.e., Radio Galaxies (RGs) an d Steep Spec- 
trum Radio Quasars (jAbdo et al. 2010at ). In spite of 
their small number, these 7-ray emitters are extremely 
appealing sources, as they offer a different perspective to 
approach the high energy phenomena. 

For example, the study of the Spectral En- 
ergy Distribution (SED) of Fan aroff-Riley I (FRI ) 
RGs su ch as NGC 1275 (lAbdo et al. 2009al) . 
M87 (lAbdo et al 2009bh and NGC 6251 
(|Migliori et al. 201 If ) has clearly evidenced the dif- 
ficulty of the one-zone homogeneous Synchrotron 
Self-Compton model (SSC) in describing the jet emis- 
sion within the AGN unified scheme. The SED modeling 
required low Lorentz factors (r < 3), implying bulk 
velocities of the RG jets much lower than typical values 
found in BL Lac objects (BLs). This is in contrast with 
the idea that blazars and MAGN are the same kind 
of source, only characterized by a different inclination 
angle of the jet (Urry fe Padovani 1995). As suggested 
by several authors (IGeorganopoulos fc Kazanas 20031 : 
iGhisellini et al. 20051: iGiannios et al. 2010 f ). this dif- 
ficulty can be resolved if the hypothesis of a one-zone 
homogeneous emitting region is relaxed and a structured 
jet is assumed. In these models, an efficient (radiative) 
feedback between different zones at different velocities 
can reconcile the apparent FRI-BL discrepancy. In the 
MAGN sample detected by the Fermi-LAT, there is a 
clear predominance of nearby RGs with FRI morphol- 
ogy. High power radio sources (i.e., FRIIs) with GeV 
emission are rare. 3C 111 is indeed the only FRII Radio 
Galaxy of the sample with a secure 7-ray association. 
The rarity of 7-ray counterparts of FRIIs could be due 



to intrinsic jet differences (jGrandi 20111 ). 

The exact localization of the 7-ray emitting region is 
however still an open problem. Fast GeV flux varia- 
tions, observed on timescales of hours to days (Abdo et 
al 2010c) have unambiguously attested to the compact 
nature (but not the localization) of the 7-ray source in 
blazars. Very recently, Agudo et al. (201 la, b) claimed 
that multiwavelength outbursts, observed in two BLs, 
i.e., OJ 287 and AO 0235 + 164, can occur in differ- 
ent zones along the jet, in the stationary core as well 
as in knots, even at larger distances (> 14 pc) from the 
black hole. Furthermore, in the Flat Spectrum Radio 
Quasar PKS 1510—089, the complexity of the radio to 
7-ray variability was interpreted as being due to differ- 
ent events occurring a long the jet as a knot propagates 
(|Marscher et al. 20101) . Constraining the black hole-7- 
ray source distance is a critical issue. Indeed it has a 
strong impact on the physical models invoked to explain 
the high energy emission (Tavecchio et al. 2011). 

In Radio Galaxies, the situation is even more compli- 
cated. Time variability is difficult to detect at GeV ener- 
gies, because of poor statistics. In addition, the discov- 
ery of 7-ray emission from the radio lobes of Centaurus 
A shows that Ge V photons can also originate from ex - 
tended structures (jAbdo et al. 2010cUlYang et al. 2011 ). 
In only one case, the FRI RG NGC 1275, flux variations 
on time scal es of months were clearly detected in the L AT 
light curves ijAbdo et al. 2010al iKataoka et al. 20101) al- 
lowing the compact nature (but not the localization) of 
the 7-ray emitting region to be asserted. Only in the case 
of the FRI Radio Galaxy M87, some clues on the loca- 
tion were provided by the very high energy (VHE) pho- 
tons, detected by the ground based Cherenkov telescopes 
(Aharonian et al. 2003, 2006; Acciari et al. 2008, Albert 
et al. 2008). The TeV flare in 2005 was simultaneous to 
a strong X-ray burst associated with the innermost knot 
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in the jet (Harris et al. 2006). In 2008, the coincidence 
of the VHE and radio flares, monitored with the Very 
Long Baseline Array (VLB A), confined the VHE photon 
production within ten Schwarzschild radii (Acciari et al. 
2009). A subsequent TeV flare in 2010 was related to 
a flux increase of the core in X-rays (Raue et al. 2011). 
The difficulty in the data interpretation does not allow an 
unequivocal spatial identification of the emitting region. 

In this paper we present the first localization of a 7- 
ray flare in a FRII Radio Galaxy. The source, 3C 111, 
was observed by Fermi in a high state in late 2008, while 
a radio-optical- X-ray campaign was ongoing (Chatterjee 
et al. 2011, hereafter CH11). The mm-to-7-ray light 
curves combined with high resolution VLBA images in- 
dicate that the GeV flare occurred during the ejection of 
a bright plasma knot from the radio core. 

2. THE 3C 111 7- RAY FLARE 

2.1. The source 

3C 111 (z=0.0491; Eracleous et al. 2004), opti- 
cally classified as a Broad Line Radio Galaxy (BLRG) 
shows a typical FRII morphology. It exhibits radio 
lobes, a strong core and a prominent 120" one-sided 
j et, terminating in the hot spot of the northeast lobe 
(|Linfield fc Perlev 1984) . At a distance of 196 Mpc, the 
linear extension of the jet is 114 kpc (1"=948 pc) in pro- 
jection, assuming H Q — 71 km s -1 Mpc -1 , Q m — 0.27, 
and J7a = 0.73. 

VLBA observations reveal a one-sided core-jet structure 
characterized by an intrinsic half o pening angle 4> ~ 3 ° 
and an inclination angle 9 ~ 18° (jJorstad et al. 2 005). 
The milliarsecond (mas) jet contains multiple super- 
lumin al components with apparent velocities f3 app = 
2 — 6 (Kadlcr et al. 20081), ejected approxima tely every 
5 months from the core (|Jorstad et al. 20 05). The ac- 
cretion disk of this source is very efficient in convert- 
ing gravitational power into radiation, as attested by the 
presence of an iron Ka line detected by several X-ray 
satellites (CH11, Ballo et al. 2011). 

3C 111 was one of the few sources in the MAGN sample 
proposed as a possible 7-ray counterpart of an EGRET 
detection (jHartman et al. 2008 and references therein ) 
and subsequentl y confirmed by Fermi. In the first LAT 
catalog of AGN lilLAC: Abdo et al. 2010bh . based on 11 
months of data, the source was associated with the low 
latitude 7-ray source 1FGL J0419. 0+3811 (detection sig- 
nificance a — 4.3) with a high probability (P = 87%). 
3C 111 was still detected with a similar significance 
after 15 and 24 mo nths of survey (|Abdo et al. 2010at 
iKataoka et al. 2011). It is not listed in the second 
LAT AGN catalog ( Ackermann et al. 201 1| ). where only 
sources with a detection > 5er are included. Indication of 
possible time variability was first suggested by Hartman 
et al. (2008). They noted that the source only occasion- 
ally became bright (-F100 > 10~ 7 phot cm -2 s -1 ) and 
detectable by EGRET, suggesting a low duty cycle for 
the 7-ray emission. The Fermi light curve, based on 15 
months of data, revealed a similar behavior. 3C 111 was 
detected with a sig nificance larger tha n 3a only in 2008 
August- November (|Abdo et al. 2010al ). 

2.2. The 24-month j-ray light curve 



Taking advantage of the recent improvement of the in- 
strument response functions (IRFs) and the background 
models, we produced a new 0.1-100 GeV light curve of 
3C 111 from 2008 August to 2010 August using a bin 
time of 2 months. The Fermi- LAT ScienceTools software 
(version v9r23pl) was used with the P7SOURCE_V6 
set of instrument response functions. The time inter- 
vals when the rocking angle of the LAT was greater 
than 52° were rejected and a cut on the zenith angle 
of 7-rays of 100° was applied. The model for which 
we calculated a binned likelihood is the combination of 
point-like and diffuse sources located within a region of 
interest (Rol) having a radius of 15° and centered on 
3C 111. For each point-like source the spectral slope 
was frozen to the best fit val ue provided by the Ferm i- 
LAT Second Source Catalog (|2FGL: Nolan et al. 2012h . 
The Rol model also includes sources falling between 15° 
and 20° from the target source, which can contribute to 
the total counts observed in the Rol due to the energy- 
dependent size of the point-spread function of the instru- 
ment. For these additional sources, normalizations were 
also fixed to the values of the 2FGL catalog. The spec- 
tral model used for 3C 111 is a power law (F=kE _r ) 
with the spectral slope fixed to the value reported in 
1LAC paper r = 2.6 ± 0.2. Note that the T uncer- 
tainties have very negligible effect on the flux calcula- 
tion. A Galactic emission model (gal_2yearp7v6_v0.fits) 
and an extragalactic and instrumental background 
(isotrop_2year_P76_source_v0.txt) were used to model 
the background. We evaluated the significance of the de- 
tection given by the Test Statistic TS=2Alog(likelihood) 
between the models with and without the source (Mat- 
tox et al. 1996). When the TS was less than 10, the 
bin flux value at F > 0.1 GeV was replaced by a 2a 
upper limit, calculated by finding the point at which 
2Alog(likelihood)=4 when increasing the flux from the 
maximum-likelihood value. 

The 24-month light curve confirms the high state of 
3C 111 in late 2008 (Fig. 1-left panel): the source ap- 
pears to spend not more than two months flaring and 
then drops below the LAT sensitivity threshold. In order 
to better constrain the 7-ray flare shape, a more dense 
LAT time sampling from 2008 August to December was 
then performed. A light curve was generated by dividing 
the time interval into bins of 1 month duration and re- 
peating the likelihood analysis for each interval (Fig. 1- 
right panel). The source is detected with a significance 
of ~ 3cr (TS=9.7) only on one occasion, indicating that 
the GeV flare peaked between October 4 and November 
4. It is probable that the major 7-ray activity was at the 
beginning of October. A likelihood analysis performed 
using a sampling of only 15 days shows that the largest 
Test Statistic value (TS=7) is observed between October 
4 and 19. 

2.3. Localization of the j-ray dissipation region 

As a multiwavelength campaign of 3C 111 (CH11) was 
ongoing in the same period monitored by LAT, it was 
possible to compare the 7-ray data with mm-optical- 
X-ray light curves (Figure 2). The simultaneity of the 
flare is impressive: the source luminosity was increasing 
from mm to X-ray frequencies exactly when the great- 
est flux of 7-ray photons occurred. This is a clear in- 
dication of the co-spatiality of the event. As attested 
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by the TeV-radio campaigns organized for the RG M87 
(Acciari et al. 2009), a fruitful approach to identify the 
zone where high energy photons are produced is to fol- 
low the jet dynamics at mas scales. We therefore in- 
vestigated the VLBA images at 43 GHz made available 
by the blazar research group at Boston UniversitjQ and 
discussed in CH11. Looking at Figure 6 of CH11, it 
appears evident that a deep modification of the core- 
jet system took place in late 2008. As deduced by the 
same authors, on October 29 (with an uncertainty of 
~ 26 days) a new component was ejected from the sta- 
tionary core. In order to quantify the jet dynamics, we 
analyzed the VLBA images at four epochs, before (2008 
September) and after (2008 November-December, 2009 
January) the appearance of the new blob. To minimize 
the free parameters, the model fitting of the (u,v) data 
were performed with DIFMAP (Shepherd 1997) using 
the simplest models (i.e., delta, circular Gaussian, el- 
liptical Gaussian). The model components were added 
to the residual images with a restored beam size of 0.1 
mas (Figure 3). This is the angular resolution in the di- 
rection of the jet propagation, obtained with a suitable 
weighting of the visibilities. The brightest component 
was identified as the stationary core, while the western- 
most component, at a distance of ~ 0.1 mas from the 
core, was associated with a possible counter-jet. The 
new knot was emerging in November and appears clearly 
separated from the core in the successive epochs. It is be- 
yond the aim of this paper to attest the real nature of the 
westernmost feature. However, we note that our analysis 
of 8 epochs, from 2008 June to 2009 February, indicates 
that the candidate counter-jet does not show any signif- 
icant systematic displacement. In the assumption that 
the counter-jet is stationary with respect to the core, we 
can estimate the position errors from its position scat- 
ter. Using the mentioned 8 epochs, the scatter is 0.005 
mas. The flux densities of the different jet components, 
when close to each other, are difficult to disentangle even 
with the model fitting procedure used here. For this rea- 
son, we preferred to follow the core evolution, considering 
the variation of the total flux density, including the con- 
tribution of the stationary feature, the counter-jet and, 
when separated, the new ejected knot. In 2009 January, 
a small component between the head of the ejected blob 
and the core is evident. Its flux was also included, given 
the difficulty in establishing its nature. It could be a new 
emerging knot or simply a disruption/elongation of the 
component ejected in 2008 (see also GroBberger et al. 
2011). The speed of the knot, considering the blob-core 
distance at three different epochs (2008 November - 2009 
January), is 1.07±0.05 mas yr _1 ((3 app = 3.3c) and the 
ejection time, obtained by extrapolating the linear (core- 
knot distance vs time) fit, is T =2008.80±0.01. Our T 
is in agreement, within the uncertainties, with the value 
of CH11, while our j3 app is smaller, because we considered 
only the 3 epochs closest to the ejection. It is interesting 
to notice that, assuming an inclination angle of 18° and 
apparent velocity (3 a pp — 3.3c, the arm length ratic0 is 

1 |http: / /www.bu.edu/blazars/l 

^ The arm Ienght ratio is the ratio between the apparent length of the 
approaching (di) and receding {(I2) jets, defined as Q— — , 
being (3 the bulk jet velocity in units of the speed of light (Longair & 
Riley 1979) 



about 20. As a consequence, the observed approaching 
components extending ~ 2 mas from the core (Figure 3) 
should be compressed to ~ 0.1 mas on the receding side. 
This latter value is consistent with the unresolved size of 
the counter-jet in the VLBA images. 

The light curve of the core is shown in Figure 4. The 
brightness increases from August to November, follow- 
ing the evolution of the multiwavelength flare (Figure 2), 
then decreases when the new radio component becomes 
visible (Figure 3). 

Finally, we note that two other knots could have been 
ejected by the core at the beginning of 2009 (CH11). 
These minor events caused no significant amplification 
of the GeV radiation. It seems that only exceptional jet 
perturbations can boost the gamma-ray flux above the 
LAT sensitivity threshold. 

3. CONCLUSIONS 

Marscher et al. (2002) in an enlightening paper on the 
RG 3C 120 noted that a dip in the X-ray light curve is 
generally associated with an expulsion of a bright super- 
luminal radio knot. It was suggested that, as in micro- 
quasars, the inner regions of an efficient accretion disk 
could be accelerated, causing a shock front streaming 
along the jet. 

3C 111 is quite similar to 3C 120. They are both 
BLRGs probably hosting a disk-corona system (Haardt 
& Maraschi 1991). CH11, exploring the radio optical 
and X-ray light curves of 3C 111 noted, also for this 
source, a strong suppression of X-rays followed by the 
appearance of a new component in the high resolution 
radio maps. They interpreted the flux increase, quasi- 
simultancously observed at mm-optical and X-ray fre- 
quencies after the dip, as the passage of a new knot 
through the jet stationary core (see also Bell & Comeau 
2011 for a different interpretation). Incidentally, we note 
that the Suzaku X-ray satellite observed 3C 111 before 
the flare, during the dip, on 2008 August 22, and XMM- 
Newton observed it after t h e blob ejection in 20 09 Febru- 
ary 15 (jBallo et al. 20111 : iTorresi et al. 2012h . As ex- 
pected, there is a difference in the 2-10 keV flux by about 
a factor 2 between the observations. 

An absorption blueshifted feature, associated with the 
resonant Fe XXVI Lya transition, was revealed by Suzaku 
but not by XMM-Newton. This absorption line could 
originate in a fast disk- wind (Ballo et al. 2011 and ref- 
erences therein). If this detection will be confirmed, the 
episodic presence of hard X-ray absorption features could 
represent a signature of an ongoing disk-jet perturbation. 

Our analysis of 24-month Fermi data reveals that the 
production of GeV photons is related to these phenom- 
ena. First of all, the short time interval of 3C 111 de- 
tectability (At ^ 30 — 60 days) by the LAT denotes the 
compactness of the emitting region. Causality arguments 
show that it cannot be larger than R < AtSc < .1 pc , 
if a Doppler factor of 5 ~ 3 (jJorstad et al. 200 5) is as- 
sumed. This immediately allows us to exclude the radio 
lobes as the main source of 7-rays. Note that assuming 
an extreme time variability of 15 d reduces the radial 
extent to 0.04 pc. In addition, 3C 111 was detected 
by the LAT exactly when the mm, optical and X-ray 
fluxes were rising. The temporal coincidence of the mm- 
optical X-ray and 7-ray flares unambiguously implies a 
co-spatiality of the event. Since the outburst of pho- 
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tons is directly connected to the ejection of a new radio 
knot, it is natural to localize the 7-ray source in the ra- 
dio core region, the size of which cannot be larger than 
0.3 pc, corresponding to the VLBA angular resolution of 
~ 0.1 mas. 

The distance between the stationary radio structure 
and the black hole (d&h— c ) is unknown. However, a 
weak VLBA feature is seen at a distance of ~ 0.1 mas 
from the core (d c _ C j), just on the opposite side of the 
ejected bright knot. If this is the counter-jet, as we sug- 
gest, the black hole necessarily lies within 0.1 mas (i.e., 
dbh-c<d-c-cj)- As a consequence, the distance of the 7- 
ray dissipation region from the central engine has to be 
of the same order of magnitude or slightly larger depend- 
ing on the position of the 7-ray source with respect to 
the counter-jet. This black hole-7-ray source distance 
(~0.1 mas corresponding to 0.3 pc), is surely a conserva- 
tive upper limit. Indeed a recent estimate of the separa- 
tion between the black hole and the base of the jet in the 
RG M87 (Hada et al. 2011) indicates shorter distances 
between these components: the stationary core seems to 
be at not more than about twenty Schwarzschild radii 
from the central engine. 
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Fig. 1. — Left panel — LAT light curve of 3C 111 covering the first 2 years of survey (from 2008 August 4 to 2010 August 4), considering 
a bin time of 2 months. The source was only detected in 2008 October-November. Right panel — One month light curve from 2008 August 
to 2009 January (lower box) and corresponding Test Statistic for each time interval (upper box). The source is detected with a TS=9.7 
(~ 3cr) only in 2008 October. 
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Fig. 2.— Millimeter (230 GHz), optical (R), X-ray (2.4-10 keV) and 7-ray (0.1-100 GeV) light curves. The original mm-optical- X-ray 
data are from CH11. A detailed study of the correlation among the different light curves is also presented in the same paper. The black 
arrow indicates the time when the radio knot was ejected by the core. 
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3C1 1 1 at 43.218 GHz 2008 Dec 21 



3C1 1 1 at 43.21 8 GHz 2009 Jen 24 
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FlG. 3. — VLBA images (43 GHz) at 4 epochs, mapping the emergence of a new knot from the radio core. The different jet components 
are shown as crosses (delta profile), circles (circular Guassian profile) and ellipses (elliptical Guassian profile). The sequence shows the 
appearance of a new component (blue filled circle) in 2008 November that is clearly separated from the core (red filled circle) in 2008 
December and 2009 January. 
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Fig. 4. — Flux Density evolution of the core region as a function of the time from 2008 September 10 to 2009 January 24. The flux 
density was estimated in a region including the core, the counter-jet and the emerging knot (see text). The arrow defines the ejection time 
as estimated by CH11. 



